KASCADE: Astrophysical results and tests of hadronic interaction models 
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KASCADE is a multi-detector setup to get redundant information on single air shower basis. The information 
is used to perform multiparameter analyses to solve the threefold problem of the reconstruction of (i) the unknown 
primary energy, (ii) the primary mass, and (iii) to quantify the characteristics of the hadronic interactions in the 
air-shower development. In this talk recent results of the KASCADE data analyses are summarized concerning 
cosmic ray anisotropy studies, determination of flux spectra for different primary mass groups, and approaches to 
test hadronic interaction models. Neither large scale anisotropics nor point sources were found in the KASCADE 
data set. The energy spectra of the light element groups result in a knee-like bending and a steepening above 
the knee. The topology of the individual knee positions shows a dependency on the primary particle. Though 
no hadronic interaction model is fully able to describe the multi-parameter data of KASCADE consistently, the 
more recent models or improved versions of older models reproduce the data better than few years ago. 



1. Introduction 

The all-particle energy spectrum of cosmic rays 
shows a distinctive feature at few PeV, known 
as the knee, where the spectral index changes 
from —2.7 to approximately —3.1 (Fig. QJ. At 
that energy direct measurements are presently 
not possible due to the low flux, but indirect mea- 
surements observing extensive air showers (EAS) 
are performed. Astrophysical scenarios like the 
change of the acceleration mechanisms at the cos- 
mic ray sources (supernova remnants, pulsars, 
etc.) or effects of the transport mechanisms in- 
side the Galaxy (diffusion with escape probabili- 
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ties) are conceivable for the origin of the knee as 
well as particle physics reasons like a new kind 
of hadronic interaction inside the atmosphere 
or during the transport through the interstellar 
medium. 

Despite EAS measurements with many experi- 
mental setups in the last five decades the origin of 
the kink is still not clear, as the disentanglement 
of the threefold problem of estimate of energy and 
mass plus the understanding of the air-shower de- 
velopment in the Earth's atmosphere remains an 
experimental challenge. For a detailed discussion 
of the subject see a recent review [T. 

The multi-detector system KASCADE (KArls- 
ruhe Shower Core and Array DEtector) ap- 
proaches this challenge by measuring as much as 
possible redundant information from each single 
air-shower event. The multi-detector arrange- 
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Figure 1. Primary cosmic ray flux and primary 
energy range covered by KASCADE and its ex- 
tension KASCADE-Grande. 

ment allows to measure the total electron and 
muon numbers of the shower separately using 
an array of shielded and unshielded detectors at 
the same place. Additionally muon densities at 
further three muon energy thresholds and the 
hadronic core of the shower by an iron sampling 
calorimeter are measured. 

In the following we present the main results 
of KASCADE, in particular the reconstruction 
of energy spectra of single primary mass groups 
and approaches of correlation analyses to test 
the hadronic interaction models and to find con- 
straints for the improvement and development of 
the next generation of the models. These tests 
provide complementary information to the data 
of present accelerator experiments, as air-shower 
data are sensitive to higher energy interactions 
and to a different (extreme forward direction) 
kinematic region. 

2. The KASCADE experiment 

The KASCADE experiment [2] measures 
showers in a primary energy range from 100 TeV 
to 80PeV and provides multi-parameter mea- 
surements on a large number of observables con- 
cerning electrons, muons at 4 energy thresholds, 
and hadrons. The main detector components of 
KASCADE are the Field Array, the Central De- 
tector, and the Muon Tracking Detector (Fig.[2J). 

The Field Array measures the total electron 
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Figure 2. The main detector components of the 
KASCADE experiment: (the 16 clusters of) Field 
Array, Muon Tracking Detector and Central De- 
tector. The location of some stations of the 
Grande array is also displayed. The lower part 
of the figure shows a sketch of the Central Detec- 
tor with its different detection set-ups. 



and muon numbers (i?^ > 240 MeV) of the 
shower separately using an array of 252 detec- 
tor stations containing shielded and unshielded 
detectors at the same place in a grid of 200 x 
200 m^. The excellent time resolution of these 
detectors allows also decent investigations of the 
arrival directions of the showers in searching large 
scale anisotropics and, if exist, cosmic ray point 
sources. 

The Muon Tracking Detector measures the in- 
cidence angles of muons {E^ > 800 MeV) relative 
to the shower arrival direction. 
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The hadronic core of the shower is measured 
by a 300 iron samphng calorimeter installed 
at the KASCADE Central Detector: Three other 
components - trigger plane (serves also as tim- 
ing facility), multiwire proportional chambers 
(MWPC), and hmited streamer tubes (LST) 
- offer additional valuable information on the 
penetrating muonic component at 490 MeV and 
2.4 GeV energy thresholds. 

The redundant information of the showers mea- 
sured by the Central Detector and the Muon 
Tracking Detector is prevailingly used for tests 
and improvements of the hadronic interaction 
models unavoidably needed for the interpretation 
of air shower data. 

3. KASCADE results 

First, some results which are obtained nearly 
independently of influences of the simulations, in 
particular studies of anisotropics of cosmic rays 
and gamma ray search will be discussed. Second, 
results of the main analyses of KASCADE are re- 
ported, i.e. the imfolding of the two-dimensional 
electron number to muon number size spectrum 
in energy spectra of five different mass groups and 
the dependence of the results on the hadronic in- 
teraction models. Finally, specific tests of the 
hadronic interaction models using data of the 
Central Detector and the Muon Tracking Detec- 
tor of KASCADE are discussed. 

3.1. Search for anisotropies and point 
sources 

Investigations of anisotropies in the arrival di- 
rections of the cosmic rays give additional infor- 
mation on the cosmic ray origin and of their prop- 
agation. Depending on the model of the origin of 
the knee one expects large-scale anisotropies on a 
scale of 10"'* to 10~^ in the energy region of the 
knee and depending on the assumed structure of 
the galactic magnetic field. For example in Fig. O 
(upper panel) the predictions from calculations of 
Candia et al. are compared with the limits of 
anisotropy given by KASCADE results j^. The 
KASCADE limits were obtained by investigations 
of the Rayleigh amplitudes and phases of the first 
harmonics. Taking into account possible nearby 
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Figure 3. Upper part: Rayleigh amplitude of 
the harmonic analyses of the KASCADE data |3] 
(limit on a 95% confidence level) compared to the- 
ory predictions Lower part: Significance dis- 
tributions for searching point sources on the sky 
map seen by the KASCADE experiment 0. 

sources of galactic cosmic rays like the Vela Su- 
pernova remnant the limits of KASCADE al- 
ready exclude particular model predictions. But 
for a complete picture the investigations have to 
be performed with air shower samples of the dif- 
ferent mass groups which need a higher statistical 
accuracy in measurements. 

Some interest for looking to point sources in the 
KASCADE data sample arises from the possibil- 
ity of unknown near-by sources, where the deflec- 
tion of the charged cosmic rays would be small 
or by sources emitting neutral particles like high- 
energy gammas or neutrons. Due to their small 
decay lengths the latter ones are of interest for 
near-by sources only. Fig. |31 (lower panel) shows 
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Figure 4. Derived limits on the gamma flux com- 
pared to limits of other experiments and theoret- 
ical spectra from [5]. 

the distribution of significances for a deviation 
of the flux from the expected background for all 
bins of the visible sky of KASCADE. Shown arc 
the distributions for the full sample of air showers 
as well as for a sample of "muon-poor" showers 
which is a sample with an enhanced number of 
candidates of 7-ray induced events. No signifi- 
cant excess was found in both samples [H]. 

With a similarly obtained sample of enhanced 
gamma candidates an analyses is performed to 
estimate limits [7] on the diffuse high-energy 
gamma ray flux. Fig. ^ shows the flux compared 
to obtained limits from other experiments and 
theoretical spectra from 8 . By increasing the 
statistics and/or improving the selection proce- 
dures the KASCADE results will reach the sensi- 
tivity to proof such theoretical calculations. 

3.2. Energy spectra of individual mass 
groups 

The content of each cell of the two-dimensional 
spectrum of electron number vs. muon number 
(Fig. ISJ is the sum of contributions from the indi- 
vidual primary elements. Hence the inverse prob- 
lem g{y) = J K{y,x)p{x)dx with y — {Ne,N^^) and 
X = {E, A) has to be solved. This problem results 
in a system of coupled Fredholm integral equa- 
tions of the form 
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where the probability j)a is a further integral with 



Two dimensional electron 
muon (A'^^' = number of muons in 
40-200m core distance) number spectrum mea- 
sured by the KASCADE array. The lines display 
the most probable values for proton and iron 
primaries obtained by CORSIKA simulations 
employing different hadronic interaction models. 

the kernel function kA ~ ta ■ ■ sa factorized 
into three parts. The quantity r a describes the 
shower fluctuations, i.e. the distribution of elec- 
tron and muon number for given primary energy 
and mass. The quantity tA describes the trigger 
efflciency of the experiment, and sa describes the 
reconstruction probabilities, i.e. the distribution 
of reconstructed Nf. and for given true num- 
bers of electrons and muons. The probabilities 
Pa are obtained by parameterizations of Monte 
Carlo simulations for flxed energies using a mod- 
erate thinning procedure as well as fully simulated 
showers as input of the detector simulations. 

The application of the unfolding procedure to 
the data is performed on basis of two different 
hadronic interaction models (QGSJETOl [H], 
SIBYLL2.1 [ini) as options embedded in COR- 
SIKA ^Jj for the reconstruction of the kernel 
functions By applying the above described 

procedures to the experimental data energy spec- 
tra are obtained as displayed in Fig. El 

Knee like features are clearly visible in the all 
particle spectrum, which is the sum of the un- 
folded single mass group spectra, as well as in 
the spectra of primary proton and helium. This 
demonstrates that the elemental composition of 
cosmic rays is dominated by light components be- 
low the knee and dominated by a heavy compo- 
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Figure 6. Result of the unfolding procedure. 
Upper part: based on QGSJETOl; Lower part: 
based on SIBYLL2.1. 



nent above the knee feature. Thus the knee fea- 
ture originates from a decreasing flux of the light 
primary particles. This observation corroborates 
results of the analysis of muon density measure- 
ments at KASCADE which are performed 
independently of the unfolding procedure. 

3.3. Inaccuracies of hadronic interaction 
models 

Comparing the unfolding results based on the 
two different hadronic interaction models, the 
model dependence when interpreting the data is 
obvious. Modeling the hadronic interactions un- 
derlies assumptions from particle physics theory 
and extrapolations resulting in large uncertain- 
ties, which are reflected by the discrepancies of 
the results presented here. The most prominent 
difference lies in the larger contribution of heavier 
primaries in case of the SIBYLL model, especially 
at high energies. To understand the differences 



between the results based on the two different 
models and to judge the validity of the models 
in general detailed investigations of the results 
are performed. In Fig. [S] the predictions of the 
Ne and N'^J' correlation for the two models are 
given in case of proton and iron primaries. It 
is remarkable that all four lines have a more or 
less parallel slope which is different from the data 
distribution. There, the knee is visible as kink 
to a flatter N^-Nj^ dependence above Njjl' « 4.2 . 
The heavier primary contribution on the results 
based on the SIBYLL model is due to predic- 
tions of a larger ratio of muon to electron num- 
ber for all primaries. Comparing the residuals 
of the unfolded two dimensional distributions for 
the different models with the initial data set we 
conclude ^1] that at lower energies the SIBYLL 
model and at higher energies the QGSJET model 
are able to describe the correlation consistently, 
but none of the present models gives a contenting 
description of the whole data set. Also a prelimi- 
nary analysis using the FLUKA (TBI code instead 
of GHEISHA 16 as low energy interaction model 
shows no conspicuous improvement of the situa- 
tion. 

Crucial parameters in the modeling of hadronic 
interaction models which can be responsible for 
these inconsistencies are the total nucleus-air 
cross-section and the parts of the inelastic and 
diffractive cross sections leading to shifts of the 
position of the shower maximum in the atmo- 
sphere, and therefore to a change of the muon 
and electron numbers as well as to their corre- 
lation on single air shower basis. The multiplic- 
ity of the pion generation at all energies at the 
hadronic interactions during the air shower de- 
velopment is also a 'semi-free' parameter in the 
air-shower modeling as accelerator data have still 
large uncertainties. 

3.4. Tests with hadronic observables 

Arbitrary changes of free parameters in the in- 
teraction models will change the correlation of all 
shower parameters. Tests using KASCADE ob- 
servables, which are measured independently of 
such used in the unfolding procedure, may give 
further constraints, e.g. by investigating corre- 
lations of the hadronic shower component with 
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electron or muon numbers. The aim is to provide 
hints for the model builder groups how the pa- 
rameters (and the theory) should be modified in 
order to describe all the data consistently. 

The applied method here is to evaluate the 
measured data relative to simulations (including 
the detector response) of proton and iron pri- 
maries. The measurements have to lie between 
these extreme values, otherwise the simulations 
cannot describe this specific observable correla- 
tion. Direct comparisons between data and simu- 
lations are not possible due to the unknown com- 
position of the primary particles generating the 
air showers. 

These kind of tests are performed for a large 
set of interaction models employed in the simu- 
lation package CORSIKA [Ilj. In ref. dH first 
results of these tests were published investigating 
the models VENUS 4.12 [THI, SIBYLL 1.6 
and QGSJET 98 0. The general conclusion was 
that QGSJET described the data best at that 
time, whereas strong hints could be given that 
SIBYLL 1.6 generates too few muons. These re- 
sult triggered improvements of the model leading 
to the newer version SIBYLL 2.1 [TOI. Later (201, 
NEXUS 2 Eg, SIBYLL 2.1 and QGSJET 01 EH 
were investigated with the result that the differ- 
ences between the models got smaller. Whereas 
QGSJET 01 and SIBYLL 2.1 can describe now 
the KASCADE hadronic observables, NEXUS 
calculations predict too little hadronic energy 
at observation level. Present investigations 
comparing the data with DPMJET 2.55 
QGSJET 01, and SIBYLL 2.1 confirm that these 
three models can describe the hadronic observ- 
ables and their correlations with the electron and 
muon component within the sensitivity of the 
KASCADE experiment (Fig. d), at least in the 
energy range below 10 PeV. 

The tests will be continued with new models, 
e.g. QGSJET II gSl, and with better statisti- 
cal accuracy at higher primary energies. Addi- 
tionally another approach will be followed: To 
change specific parameters in one certain model 
(QGSJET 01) and to investigate the change of the 
correlations of measurable observables and 
compare them with the data. 

A different window to the hadron-muon corre- 
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Figure 7. Correlation between the reconstructed 
hadronic energy sum and the number of muons. 

lation is the investigation of the KASCADE trig- 
ger rates [27]. Here, in particular primaries of 
lower energies and their behavior in the atmo- 
sphere is approached. By this, indications for 
an underestimate of the non-diffractive inelastic 
cross-section in the models are given. 

3.5. Tests with muon densities 

The tests described in last sections are sensi- 
tive to the high-energy interaction models used 
in EAS simulations. Measurements of muon den- 
sities for different energy thresholds, here 2.4 GeV 
and 490 MeV detected by the KASCADE Central 
Detector are more sensitive to low-energy interac- 
tion models like GHEISHA [HI, UrQMD 1.1 EHI, 
or FLUKA ^ (see also ref. EH). The ratio of 
muon densities measured at fixed core distances 
are determined by the muon energy spectrum 
in air-showers, which is a different approach to 
test the models than investigating the total num- 
ber of muons 29 . Fig. |S1 displays the ratio of 
the muon densities for two thresholds in depen- 
dence of the muon number (truncated muon num- 
ber) measured by the KASCADE array {E^'' = 
240 MeV) . The range of N^^ corresponds to the 
primary energy range of 10^ — 10^ GeV. Compared 
is the distribution of the data with predictions 
(including full detector response and reconstruc- 
tion procedures) for various high-energy and low 
energy model combinations. None of the mod- 
els can reproduce fully the measurements, espe- 
cially at higher muon number, i.e. higher pri- 
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Figure 8. Distributions of the muon density ratio 
for KASCADE data in comparison with predic- 
tions for various model combinations. 

mary energy, but the behavior of the UrQMD 
model seems to be more consistent as FLUKA 
or GHEISHA. Next generation of CORSIKA wiU 
include QGSJET II [23] as a new model, which 
show in first calculations a significantly different 
behavior of the muon component. 

3.6. Tests with muon pseudorapidities 

Another approach to test the hadronic interac- 
tion models via the muon component is possi- 
ble due to the excellent angular resolution of the 
KASCADE Muon Tracking Detector (w0.35°). 
Measuring the relative angles r and p between 
single shower muons and the shower axis the pseu- 
dorapidity rj — In-^ k, —ln(y^T^ + can be 

calculated. Fig. |^ shows the correlation of the 
pseudorapidity of the muons measurable in the 
MTD and of the parent hadrons [30]. The ra- 
pidity distribution of hadrons generated in high- 
energy hadronic interactions is still an open ques- 
tion and an important parameter for the model 
building. 

3.7. Tests of the shower development 

Investigating further shower observables at 
KASCADE like the muon arrival time distribu- 
tions jSJjj the electron and muon lateral distribu- 
tions or the muon production heights |321 enables 
more to scrutinize the shower development itself 
rather than the hadronic interactions. In sum- 
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Figure 9. Pseudorapidity distribution predicted 
by CORSIKA/QGSJET simulations for muons 
detectable in KASCADE and for the parent 
hadrons of the same muons. 

mary, within the given sensitivity of KASCADE 
CORSIKA describes the shower development rea- 
sonably well and no significant deviations were 
found. 

4. Conclusions and outlook 

Measurements with the multi-detector setup 
KASCADE provide plenty of high quality data to 
investigate the physics of the knee in the cosmic 
ray energy spectrum. By analyzing these data 
no significant large scale anisotropy is found in 
the energy region of the knee and no point source 
candidates could be identified, neither for charged 
cosmic rays nor for a sample with enhanced 
number of gamma-ray candidates. Concerning 
energy reconstruction, conclusive evidence has 
been given that the knee is caused by a decrease 
of the flux of light primaries, where the knee po- 
sitions show a dependence on the primary mass 
group. Systematic uncertainties for the estimate 
of the elemental composition are dominated by 
the inadequacy of the hadronic interaction mod- 
els underlying the reconstruction of energy spec- 
tra of single mass groups. Hence, still there are 
only weak constraints for detailed astrophysical 
models to explain the knee in the primary cosmic 
ray energy spectrum. 

Indications for the inadequate description of 
the hadronic interactions at the atmosphere are 
also given by additional KASCADE data anal- 
yses taking the advantage of the multi-detector 
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information, i.e. investigations of the hadron 
component in air-showers or of muon properties 
measured for different muon energy thresholds. 
These investigations of observable correlations 
have shown that none of the present hadronic 
interaction models is able to describe all the 
KASCADE data consistently (on a level of a few 
percent). Recently some efforts are made to sam- 
ple the information from accelerator experiments 
and cosmic ray investigations to improve the 
hadronic interaction models. 

The multi-detector concept of the KASCADE 
experiment has been translated to its extension 
KASCADE-Grande 34 , accessing higher pri- 
mary energies up to 10^* eV which will prove the 
existence of a knee-like structure for heavy el- 
ements. In future, by having the data of the 
KASCADE-Grande experiment and by further 
improving the hadronic interaction models better 
constraints especially at higher primary energies 
are expected. Thus cosmic ray physics at energies 
around the knee remains a vital field of research 
with high scientific interest. 
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